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We present an analytic model to calculate the atomic scattering length near a Feshbach resonance 
from data on the molecular binding energy. Our approach considers finite-range square-well poten- 
tials and can be applied near broad, narrow, or even overlapping Feshbach resonances. We test our 
model on Cs2 Feshbach molecules. We measure the binding energy using magnetic-field modulation 
spectroscopy in a range where one broad and two narrow Feshbach resonances overlap. From the 
data we accurately determine the Cs atomic scattering length and the positions and widths of two 
particular resonances. 
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I. INTRODUCTION 

In experiments on ultracold quantum gases, control of 
the atomic scattering length near Feshbach resonances 
has become apowerful tool to explore different interac- 
tion regimes [l| . Prominent examples include the implo- 
sion and explosion of Bose- Einstein condensates 0] , the 
creation of ultracold molecules 0, [|| , strongly interact- 
ing Fermi gasesjH, [f| , and the observation of three-body 
Efimov states Moreover, Feshbach resonances have 
been proposed as a very sensitive probe to the variation 
of fundamental constants Q. 

For all such experiments precise knowledge of the 
atomic scattering length a is desirable. Usually a multi- 
channel molecular potential model is employed to calcu- 
late the scattering length. This approach requires exper- 
imental input, as usually provided by photoassociative or 
Feshbach spectroscopic data [§, QjJ EH • It captures the 
global scattering behavior of the system and has a predic- 
tive power on Feshbach resonance positions with an un- 
certainty of typically a few hundred milligauss. In many 
experimental situations, however, precise knowledge of a 
is needed in the vicinity of specific Feshbach resonances 
at a level not reachable with this global approach. 

We here consider the situation in which the scatter- 
ing length a near a specific resonance is calculated us- 
ing molecular binding energy data. Nowadays power- 
ful experimental methods are available to measure the 
binding energy of Feshbach molecules with very high 
accuracy, such as radio- frequency and microwave spec- 
troscopy [H, [H, G3, EH, EH- For the relation between 
the binding energy E\> and a a simple analytic form can 
be given in the universal limit of very large positive val- 
ues of a. It is given by E\> = h 2 /(2m r a 2 ), with m r be- 
ing the reduced mass. In general, the relation between 
molecular binding energy and scattering length is more 



complex, depending also on short-range parameters of 
the molecular potentials. In principle one may also use a 
multi-channel calculation to precisely derive a from mea- 
surements of Ey>. When the specific scattering properties 
arc needed near a particular Feshbach resonance, a much 
more simple model can be applied. 

In this Article, we describe a versatile analytic model 
to directly convert molecular binding energy data into 
atomic scattering lengths near Feshbach resonances. Our 
approach considers square- well model potentials and can 
be applied to any broad or narrow Feshbach resonance, 
or even to the case of overlapping resonances. We demon- 
strate our model on Cs2 Feshbach molecules. Cesium is 
an excellent candidate to test scattering models because 
of its rich interaction properties [l?], E3 ■ From fitting the 
binding energy data, we show that the magnetic-field de- 
pendent scattering length can be determined with high 
precision. We compare our results with the full multi- 
channel numerical calculation of Rcf. 1171. 



II. FESHBACH RESONANCE MODEL 

We employ a two-channel square well potential 
V(R) to describe interacting atoms and weakly-bound 
molecules near a Feshbach resonance. The well size is 
chosen to account for the long range behavior of the 
molecular potential, which is dominated by the van der 
Waals (vdW) interaction V{R) = -C 6 /R 6 , where C 6 is 
the vdW coefficient and R is the atomic separation. The 
associated vdW length scale and vdW energy scales are 
RvdW = ^(2m r C 6 /h 2 ) 1 ^ and E vdw = h 2 /2m r Rl dW . 

To fully capture the threshold behavior of atomic scat- 
tering with vdW interaction, the well size is chosen to be 



the mean scattering length a = 47rr(l/4) 2 i? v dw |20j |. 
where r(x) is the gamma function. This choice yields 
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FIG. 1: (color online) Two-channel square well model of Fesh- 
bach resonances. Scattering states are initially prepared in 
the entrance channel |e); the closed channel |c) supports a 
bare bound state at energy E c . Feshbach coupling is mod- 
eled by the mixing of the two channels, which is resonantly 
enhanced when the bound state is near the scattering state. 
The interaction range is set to be the mean-scattering length 
of the atomic van der Waals potential, see text. 



the same binding energy of weakly-bound molecules in 
the vdW potential, namely, Ef, = h 2 /2m r (a — a) 2 for 

Eb <C -Evdw [13] ■ 

In our model, scattering atoms are initially prepared 
in one spin configuration, called the "entrance channel" 
\e). A second "closed channel" |c) supports a molec- 
ular bound state, see Fig.[TJ The quantum state of 
an atomic pair with energy E is described as \ip) = 
ip c (R)\c) + tp e (R)\e), which satisfies Schrodinger's equa- 
tion: 
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where, using a matrix representation, we write 
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Here for R < a, we assume the attractive potential 
can support multiple molecular states, that is, V e ,V c 3> 
.EvdW, and fifl induces Feshbach coupling between the 
channels. For R > a, entrance and closed channel thresh- 
olds are set to be E = and E = oo, respectively. See 
Fig.ffl 



A. Scattering state above threshold and scattering 
length 

For E = h 2 k 2 /2m r > 0, the solution of Eq. (JTJ) corre- 
sponds to an eigenstate lying above the entrance-channel 
threshold and is given by 
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where A is a constant, h 2 q 2 _/2m r = E 
\{V e -V c )$ec26, tan 29 = 2Ml/(V e -V c ] 



for R < a (2) 
(3) 

\{V e + V c ) ± 



+) = cosO\e) + 
The scattering 
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smt)\c) and | — ) = —single) +cos6*|c). 
phase shift r\ can be derived from the continuity condition 
of the wavefunction at R — a, which gives 
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(4) 

tan (ka + rf) tan (q+a) tan(c7_a) 

Potential depths V c and V e can be chosen such that 
in the weak coupling limit # ^ 0, the scattering phase 
shift approaches the off-resonant value rjbg and the closed 
channel bound state is at E c . These two conditions, given 
by 



(5) 
(6) 



tan(g + a) tan (ka + rybg) 



sin \ h~ 1 yj2m r (E c + V c )a) =0 



allow us to eliminate q+ and q- in Eq. Assuming 
further that only one bare bound state is near the con- 
tinuum \E C \ <C \/V c -E v dw, and the Feshbach coupling is 
weak \9\ « HSl/\V e - V c \ < 1, we get 



cot(fca + rj) = cot(ka + ?7bg) — 



r/2 

kaE r 



(7) 



where we have introduced the Feshbach coupling strength 
T/2 = 26 2 V C . 

Scattering lengths can be derived from the threshold 
relation a — — lhrik^o(k cot rf^ 1 . From Eq. ([7]), we get 



1 



1 



r/2 

aE c ' 



(8) 



a — a abg — a 

where <Jb g = — limft_»o(^ cot rj^g)^ 1 is the background 
scattering length. 

Magnetic tunability of Feshbach resonances can be 
modeled by a linear Zeeman shift of the bare bound 
state as E c = Sfi(B — B c ), where B is the magnetic field, 
Sfi is the relative magnetic moment of the two channels, 
and B c is magnetic field that tunes the bare state to the 
entrance channel threshold. This magnetic field depen- 
dence allows us to rewrite the scattering length in the 
standard resonance form a = ab g [l — A/ (B — £>o)]- Here 
the resonance width A and position Bq are given by 



A = Sfj, 1 a-, 
B = B C -PA, 

where a = (a^ g — a) 2 /ab g a and = 
dimcnsionlcss parameters. 



(9) 
(10) 

a b g / (a-bg — a) are 
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B. Bound molecular state below threshold 

We assume that a bound eigenstate \tp m ) is located 
near and below the continuum with binding energy Eb = 
ti 2 k 2 a /2m r . Solving Eq. ([T]), we obtain the wave function 
as 



\tp m {R > a)} oc \e), 

sin(q + R) A m sm(q-R) 
\ip m {R<a)) oc 1+) + - 



R 



R 



where A rn is a constant and q± = (q\ — k 2 — fc 2 ^) 1 / 2 . 

The binding energy Eb can be determined from the 
boundary conditions at R = a. Following the same 
calculation as that leading to Eq. J5]), we find that 
for large background scattering length |ab g | 3> a, the 
wave-number k m determining the binding energy Eb = 
h 2 k 2 j2m r satisfies 
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(11) 



In the limit when the Feshbach coupling approaches 
zero r — ► 0, the eigenenergies reduce to Eb = —E c and 
Eb = h 2 /2m r (abg — a) 2 (for ab g > a), which correspond 
to the uncoupled bare states in the closed and entrance 
channels, respectively. 

One can now convert the binding energy into the scat- 
tering length with the following strategy. First, the bind- 
ing energy data are fitted with Eq. (jTTJ) . The fitting pa- 
rameters are Sfi, T, ab g , and B c . Usually the mean scat- 
tering length a is fixed because of sufficient knowledge of 
Cq. In some cases, also ab g or are known so that they 
can be kept fixed during the fit procedure. Second, one 
determines the width A and the location of the Feshbach 
resonance Bq from the fitting parameters using Eq. (J9]) 
and (fTUjl . 



C. Multiple overlapping Feshbach resonances 

The above discussed model can be generalized to the 
case of multiple, overlapping Feshbach resonances, which 
commonly occurs in cold collisions of alkali atoms and 
chromium atoms. 

Let us assume that TV weakly-bound states are sup- 
ported by N independent closed channels. If the mixing 
between the closed channels and the entrance channel are 
weak and |ctbg| 3> o>, we can extend Eq. ([T]) to model the 
scattering length and molecular energy. The results arc 
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where Ei is the energy of the ith bare bound state and Ti 
is its Feshbach coupling strength to the entrance channel. 
The results presented here bear similarity to those from 
the classic Fano theory on the coupling of discrete states 
to a continuum (l8j . 

Assuming the bare states can be linearly tuned mag- 
netically, namely, Ei = Sfii(B — B c ^), where S^i is the 
relative magnetic moment and B c ^ is the crossing field 
value of the i-th bare bound state, we can express the 
scattering length as a/ab g = 1 — P(B)/Q(B), where P(x) 
and Q(x) are polynomials of {N — 1)- and iV-th order, 
respectively. If we further assume all B c /s are distinct 
and all /li's have the same sign, the scattering length 
expression can be factorized as: 
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(14) 



Here we define Bq^ as the i-th lowest pole of a/ab g , and 
B* the i-th lowest zero. This form is more convenient for 
experimental research since zeros and poles of the scatter- 
ing length can be identified experimentally. Accordingly, 
the width of the i-th Feshbach resonance can be defined 
as Aj = B* - B 0>i . 



III. BINDING ENERGY MEASUREMENTS 
AND DETERMINATION OF THE SCATTERING 
LENGTH 

A. Measurements of the binding energy 

We measure the binding energies of ultracold Cs2 Fesh- 
bach molecules in a magnetic field range where a weakly 
bound s-wave molecular state overlaps with a d-wave and 
a g-wave state [l9|. The latter two states are connected 
to the d-wave and <?-wave Feshbach resonances located at 
about 48 G and 53.5 G, respectively [l7[. The molecu- 
lar spectroscopy is performed by using a weak oscillating 
magnetic field. This technique, previously applied to a 
85 Rb atomic sample [l5| and to a 85 Rb- 87 Rb mixture 
[ill, has the advantage to determine binding energies of 
weakly bound molecules to within a few percent. The 
oscillating magnetic field stimulates a transition between 
an atom pair and a molecule for a modulation frequency 
resonant to Eb/h [IHH^, where h is the Planck constant. 

Our molecular spectroscopy is performed using ultra- 
cold Cs atoms in an optical trap. A detailed description 
of our experimental procedure can be found in Ref. [23| . 
After several pre-cooling stages, the atoms are fully po- 
larized in their absolute ground state \F = 3, mp = 3). 
We evaporatively cool the atoms in the optical trap at 
a magnetic field of 27 G [24| , where the two-body scat- 
tering length is large enough (a ks 450ao ; with ao the 
Bohr radius) to provide fast thermalization. During the 
evaporation sequence, an extra magnetic field gradient is 
applied to levitate the atoms against the gravity [13, H^] . 
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FIG. 2: (color online). Example of a loss resonance at 48.2 
G. The losses are due to the coupling between atom pairs and 
d-wave molecules. The solid line is a gaussian fit to the data. 
We find the resonant frequency to be v — 93.0(3) kHz and the 
width 5.5(6) kHz. The inset shows the typical magnetic field 
sequence used to associate Cs2 molecules with an oscillating 
magnetic field. 
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FIG. 3: (color online). Binding energy of cesium molecules 
near three Feshbach resonances as a function of the magnetic 
field. Zero energy corresponds to two Cs atoms in the ab- 
solute hyperfine ground-state sublevel \F = 3,rriF =3). The 
measurements are shown as open circles. The fit (solid line) 
is based on Eq. (|13|l . see text. The inset shows an expanded 
view in the region of the two d- and g-wave narrow resonances. 
The error bars refer to the statistical uncertainties. 



We typically obtain a ncarly-dcgcncratc sample of 8000 
Cs atoms at T = 100 nK. Typical values for the density 
and phase space density arc 1.2 x 10 12 cm" 3 and 0.2, 
respectively. 

We measure the molecular binding energy by mag- 
netic field modulation [l5|, [2^. The inset of Fig. [5] shows 
the magnetic-field sequence. After ramping the mag- 
netic field within 3 ms to a desired value in the range 
of interest, we add for typically 1 s a small sinusoidal 
magnetic field B(t) = Asin(27wi), where the modula- 
tion amplitude A typically ranges from 100 to 600 mG. 
We then ramp back the magnetic field to initial value and 
wait 100 ms, during which the associated molecules de- 
cay by inelastic collisions with the atoms. The magnetic 
field is then rapidly turned off to image the atoms with 
standard absorption imaging technique. If the modula- 
tion frequency v is close to the binding energy value, 
Eh ~ hv, the oscillating magnetic field couples atom 
pairs to molecules [L'J . The atoms are converted into 
molecules, and the resulting association is detected as 
resonant loss in the atom number. A typical loss signal 
is shown in Fig.[5J The resonant frequency is determined 
to within a few percent by fitting the data with a Gaus- 
sian profile [29[ . 

We determine the resonant loss frequency for different 
values of the magnetic field. Figure[3]shows the measured 
molecular spectrum (open circles) in the magnetic field 
range between 25 G to 60 G. Here, three weakly bound 
molecular states, namely s-, d-, and <?-wave bound states, 
overlap near threshold [HJ, [2l| . We usually observe losses 
of Cs atoms up to 50%. The largest loss rate and hence 
conversion efficiency is observed when g-wave molecules 
are associated. This observation is consistent with the 
results of Ref. [13], which predict higher conversion effi- 
ciencies for larger relative magnetic moments 6[A. 



B. Determination of the scattering length 



We now convert the binding energy data into scatter- 
ing length. For Cs atoms, the vdW energy is -E v dw,Cs = 
h x 2.7 MHz, and the mean scattering length a = 95.7<2b- 
We apply our model in the regime Et, < £vdw,Cs by 
fitting the binding energy data with Eq. (fl"3|) for over- 
lapping resonances. We consider three resonance terms, 
which correspond to couplings to the s-, d- 1 and g-wave 
molecular states. Each term contains three parameters: 
Feshbach coupling strength Ti, bare state crossing posi- 
tion B Ct i and the relative magnetic moment S/ii [30j ]. All 
the parameters are evaluated by the fit. We assign the 
uncertainties of the parameters by a resampling method 
[3lj . Second, we use the fitting parameters to determine 
a in the form of Eq. (|12j) . Table U summarizes the main fit 
results for the three overlapping resonances. We derive 
the zeroes and poles of the scattering length from Eq. (fi"2]) 
and (fT4|) . Our results represent the most accurate knowl- 
edge of the positions and widths of the two experimental 
important resonances near 48 and 53 G [25|, [26[ 

Figure [4] shows the magnetic-field dependent scatter- 
ing length as derived from our model and the binding 
energy data. We observe a remarkable overall agreement 
between these results and a previous multi-channel calcu- 
lation based on the knowledge of the positions of various 
Feshbach resonances [T3]. Our present method proves 
particularly powerful for the scattering length near the 
two narrow resonances. Here the previous calculations 
[l7| suffer from the large uncertainties of the positions of 
the resonances, which can even exceed their widths. 
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TABLE I: Fitting parameters for the s-, d- and p-wave Fesh- 
bach resonances, determining the scattering length in the 
magnetic field range of interest; see Fig. [3] The background 
scattering length db g = 1875<zo, the mean scattering length 
of cesium a = 95.7 as and the bare s-wave state magnetic 
moment 5/j,i = 2.50/is are set constant. Poles Bo,i and 
zeros B* of the scattering length are derived, see text. Un- 
certainties in the parentheses are statistical. The systematic 
uncertainty of the magnetic field is 10 mG. 
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r<//i(MHz) 


Sm/fiB B Ctl (G) B 0:l (G) B* (G) 


s-wv. 
d-wv. 
g-wv. 


11.6(3) 
0.065(3) 
0.0042(6) 


2.50 19.7(2) -11.1(6) 18.1(6) 
1.15(2) 47.962(5) 47.78(1) 47.944(5) 
1.5(1) 53.458(3) 53.449(3) 53.457(3) 
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FIG. 4: (color online) Scattering length of \F = 3,mF = 3) ce- 
sium atoms in the magnetic field range where three Feshbach 
resonances overlap. The solid curve shows the result of this 
work while the dashed curve represents the prediction from a 
previous multi-channel calculation fl7l |. 



IV. CONCLUSION 



In summary, we have described a new method to pre- 
cisely determine scattering lengths from binding energy 
measurements of weakly bound molecular states. We em- 
ploy a simple interaction model to extract the relation- 
ship between the scattering length and the binding en- 
ergy. Using the cesium molecule binding energy data in 
the range where three Feshbach resonances are present, 
we show that our model provides an excellent fit to the 
molecular energy structure, and allows us to precisely de- 
termine the scattering length near the resonances. The 
precise knowledge of binding energy and scattering length 
near Feshbach resonances provides an ideal starting point 
to explore few- and many-body physics by controlling 
the two-body interactions. Our model is also general 
in the sense that it can apply to any homonuclear or 
hctcronuclear atomic systems near Feshbach resonances, 
once molecular binding energy data are available. 
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